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1 Introduction 

The MARS scanner[T] is designed for the x-ray spectroscopic study of samples 
with the aid of computer tomography methods. Gomputer tomography allows the 
reconstruction of slices of an investigated sample using a set of shadow projections 
obtained for different angles. Projections in the MARS scanner are produced us¬ 
ing a cone x-ray beam geometry. Gorrect reconstruction in this scheme requires 
precise knowledge of several geometrical parameters of a tomograph, such as dis¬ 
placement of a rotation axis, x-ray source position with respect to a camera, and 
camera inclinations. Use of inaccurate parameters leads to a poor sample recon¬ 
struction. Non-ideal positioning of camera, x-ray source and cylindrical rotating 
frame (gantry) itself on which these parts are located, leads to the need for to¬ 
mograph alignment. In this note we describe the alignment procedure that was 
used to get different geometrical corrections for the reconstruction. Also, several 
different estimations of the final spatial resolution for reconstructed images are 
presented. 

2 The MARS scanner 

The MARS (Medipix All Resolution System) scanner (Figure obtained by 
the Laboratory of Nuclear Problems of the Joint Institute for Nuclear Research, 
has a camera with two detectors based on gallium arsenide (compensated with 
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chromium) sensors|2] and equipped with the Medipix3.1|3j electronics. Each de¬ 
tector has 256 x 256 pixels of 55 microns and has dimensions of 14 x 14 mm. 
The MedipixS.l electronics can simultaneously count the number of photons with 
energies above the threshold at each pixel for two different thresholds. Gallium 
arsenide sensors have higher detection efficiency of photons at high energies in 
comparison with silicon sensors. The advantages of the scheme with single pho¬ 
ton counting include the absence of dark currents and radiation dose reduction. 
Moreover, by changing user-settable parameters in the MedipixS.l electronics, it 
is possible to select different ranges of photon energies for spectrometry. 
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Figure 1; The MARS scanner in the Laboratory of Nuclear Problems of the Joint 
Institute for Nuclear Research. 

Since the MARS tomograph was designed for biomedical research, the gantry 
with the scanning equipment (x-ray source and camera) attached to it is rotated 
around fixed scanned sample (Figure]^. The gantry rotation axis is horizontal. 
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A test sample (up to 100 mm in diameter and 300 mm length) is placed in the 
center and can be moved along the rotation axis. The camera and the source are 
placed on the gantry from the opposite sides and are able to move along the axis 
connecting each other, approaching or moving away from the sample. In order 
to enable scanning of large samples camera can also be moved in the direction 
transverse to the rotation axis and the axis connecting the source and the camera. 



Figure 2: Photo of the MARS equipment with a scanned sample. 

2.1 The MARS tomograph reference frame 

To describe the geometry of the tomograph it is necessary to dehne the reference 
frames associated with the moving parts of the tomograph (Figure]^: 

Laboratory reference frame 

The Cartesian reference frame {X, Y, Z), associated with the tomograph 
support structure. The axes are as it is shown in the hgure. 
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Sample reference frame 

The Cartesian reference frame (X^, Yg, Zg), associated with the sample, 
which can be moved along the Z axis of the laboratory reference frame. 

Gantry reference frame 

The Cartesian reference frame (Xg, Yg, Zg), associated with the gantry, and 
differs from the laboratory frame by the rotation angle 0 around the Z axis. 

Camera reference frame 

The Cartesian reference frame Xc, Yc, Zc), associated with the camera that 
can be moved along the Yg axis of the gantry reference frame. This reference 
frame is dehned by the pixel detector plane. The Y^. axis is parallel to the 
detector rows of pixels, the Zc axis is parallel to the detector columns of 
pixels. 



Figure 3: Reference frames associated with the moving parts of the tomograph 
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3 Mechanical precision of the movement of the to¬ 
mograph parts 

3.1 Measurements using an electronic feeler gauge 

Prior to the alignment studies, the mechanical precision of the movements of the 
tomograph parts has been investigated. The main goal was to measure the dis¬ 
placement of the gantry rotation axis during its rotation and to measure the dis¬ 
placement of the gantry itself along the rotation axis. To perform these measure¬ 
ments an electronic feeler gauge was used. In the first case, the positions of the 
upper points of the left and right gantry supports (red arrows in Figure left 
side) have been measured for different rotation angles. The mean value of these 
two measurements (left and right) gives an estimation of the displacement of the 
rotation axis in the local area of a measured sample. In the second case, position 
measurements have been performed from the near side and the far side of the edge 
of the left gantry support (red arrows in Figure [fright side). In this case the mean 
value of the two measurements describes gantry movement along its rotation axis. 



A + 

Figure 4; Measurement of the rotation axis Z displacement along vertical Y-axis 
(left) and measurement of the gantry movement along the Z-axis (right). 


The results of the described mechanical measurements are presented in Fig¬ 
ure 1^ and Figure as a function of the gantry rotation angle. The results show 
that in the area of a scanned object the displacement of the gantry rotation axis 
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is withing 30 microns (green line in Figure]^, while the gantry itself moves along 
Z-axis by 100 microns (green line in Figure]^. These displacement values are not 
negligible with respect to the pixel size of the detector and will lead to additional 
smearing of a reconstructed image. The results of these measurements may be 
used in the future for correcting projections of scan depending on the angle of 
rotation. 



right support left support 


mean value 


Figure 5: Relative position of upper points of the left (red) and right (blue) gantry 
supports as a function of the gantry rotation angle. The mean value (green) 
characterizes the displacement of the gentry rotation axis in the scanning area 
during rotation. 
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mean value 


Figure 6: Position measurements of the "near" (red) and "far" (blue) sides of 
the left gantry support for different rotation angles. The mean value (green) 
characterizes the displacement of gantry along its rotation axis. 

3.2 Measurements using angular encoder system 

Rotation of the gantry was studied using the optical angular encoder system REN- 
ISHAW RESR|1|. The system we used consists of the RESR ring (with the di¬ 
ameter of 200 mm) with graduations marked on it every 20 pm and the optical 
readhead with a resolution of 5 pm. The RESR ring is mounted on the shaft of 
the gantry and the readhead is attached to the tomograph frame. The angular 
resolution of this system is about 10 arc seconds. 

In the normal scanning procedure the tomograph equipment rotates around a 
hxed sample making stops to get the x-ray projection images at different rotation 
angles. The gantry with its equipment attached has a considerable moment of 
inertia, which leads to vibrations of tomograph parts at the moment it stops. The 
behavior of the tomograph at this moment was investigated with the aid of RESR 
system (Figure [^. As shown in the hgure, significant initial vibrations of the 
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system are observed. To get projections of good quality, there should be a time 
delay between the stop and the exposure of about half of a second. 



Figure 7; Angular measurements after stopping the tomograph rotation. 

The next step of the investigation is to study the repeatability of angular 
measurements at the same positions of the step motor which rotates the gantry. 
During a scan the tomograph makes stops several times in the same angular po¬ 
sitions to get projections. For successful stitching of these projections, a high 
accuracy in the repeatability of angular positions is needed. Figure shows the 
results of numerous angular measurements when tomograph is moved to positions 
0° and 352°. The results show that the design of the tomograph does provide high 
reproducibility of angular positions. 
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Figure 8: Repeatability of angular measurements when tomograph is moved to 
positions 0° and 352°. 

The gantry of the MARS scanner is driven by a step motor via an automo¬ 
tive timing belt. The weight distribution of the equipment attached to the gantry 
is non-uniform and can lead to different tension of the belt for different rotation 
positions. The following measurements have been performed to check whether the 
tomograph angular positions correspond to actual angular positions. In Figure 
the difference between measured and expected values of angular positions is pre¬ 
sented as a function of gantry rotation angle. Deviation of the measured angle is 
significant and its maximum approximately corresponds to the value of 500 mi¬ 
crons on the surface of the RESR ring. In principle, encoder measurements are 
affected not only by the ring rotation but also by the movement of the ring axis 
during its rotation. To account for this effect, independent measurements of the 
ring axis movement were performed depending on the rotation angle. The distance 
to the surface of the RESR ring was measured using SAFIBRA optical probe. The 
amplitude of this measurement does not exceed 80 microns and can not explain 
the value of the difference in measured rotation angle. The difference between the 
measured and expected values of the rotation angle can be taken into account to 
correct gantry positions and then used in the image reconstruction procedure with 
nonlinear rotation steps. 
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Figure 9: The difference between measured and expected values of angular posi¬ 
tions as a function of gantry rotation angle. 


4 Geometrical parameters requiring calibration 

In case of ideal geometry, the Zg axis of the MARS tomograph, along which the 
sample is moving, is parallel to the gantry rotation axis Zg^ parallel to the columns 
of the camera Zc axis and parallel to the Z axis of the laboratory reference frame. 
The Yc and Z^. axes corresponding to the camera pixel rows and pixel columns are 
parallel to the Yg and Zg axes, respectively. The axis along which the x-ray source 
and the camera are moving in the direction toward the sample or from the sample, 
coincides with the Xg axis. The focus of the x-ray source and the center of the 
camera are also lying on this axis, and it is crossing with the gantry rotation axis. 
The distances from this crossing point to the source SOD, the crossing point to 
the camera ODD, and from the source to the camera SDD are assumed to be well 
determined. The Yc axis Yc, along which the camera moves is parallel to the Yg 
axis. 

In reality, the geometry has some misalignment which can be expressed in 
small corrections on shifts and rotations of the tomograph parts. Some of these 
geometrical corrections can be taken into account in the reconstruction of three- 
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dimensional image of a sample, so it is necessary to know their values accurately. 
One of these parameters is the camera inclination angle with respect to the rotation 
axis, i.e. the Otiit angle between the Yc. axis and the Yg axis. We also need to 
take into account the position of the projection of the gantry rotation axis on the 
plane of the camera. It is also important to know the position of the projection 
of the source on the camera plane along the Yj, and axes. Finally, one should 
accurately measure the distance from the source to the rotation axis SOD and 
from the rotation axis to the camera ODD. 

This work is devoted to finding the camera inclination angle Otut, the de¬ 
termination of the position of the source projection on the camera plane, and the 
refinement of the distance between the source and the camera SDD. The position 
of the projection of the rotation axis on the camera plane is computed by iterative 
methods with good precision. 

5 The procedure for the MARS scanner alignment 

5.1 Phantom description 

To study the geometry of the tomograph a very simple phantom was designed. 
This is a plastic hollow cylinder of 35 mm diameter with two tungsten 50 /rm 
wires glued on the surface of the cylinder (Figure [l0|. The tungsten wires form a 
"cross" with one wire oriented along the Z-axis and the other in the perpendicular 
direction. 
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Figure 10; Phantom for the tomograph geometry calibration 


On the camera projections, the position of the cross can be determined using 
a Gaussian approximation of the shadow in the local vicinity of a crossing point. 
An example of these measurements is presented in Figure The position of 


the shadow cross was determined for the holder position -21 mm, and the gantry 
rotation angle 40° in the first case (left side) and also for the gantry angle 220° in 
the second case (right side). 
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Figure 11; Determination of the position of the shadow cross, using a Gaussian 
approximation in the local vicinity of a crossing point. 


The above mentioned cross of wires, glued to the surface of the cylinder, 
was mainly used in the procedure for determining the source position. To hnd the 
camera inclination angle, another wire was stretched in the cross section of the 


phantom, which is shown as a vertical dashed line in Figure 10 


5.2 Determination of the camera inclination angle Otiit- 

To determine the camera inclination angle atut in the YgZg plane, the projection 
of the wire stretched in the XgYg plane is used. For this purpose the shadow 
inclination angle is analyzed for the vertical wire. Because the wire itself is not 
perfectly vertical, two projections are used: direct projection, and one obtained by 
rotating the gantry 180 °. The average angle of the two positions gives the desired 
camera inclination angle in the plane YgZg (Figure [I^. 
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Figure 12: Determination of the camera inclination angle Otiit with respect to the 
gantry rotation axis in the plane YgZg plane. 

5.3 Determination of the source position along the Z(. axis 

To determine the source position with respect to the camera along the axis, we 
used the shadow cross images for the two gantry rotation angles in six Z-positions 
of the sample holder (Figure [T^. The values of rotation angles of 40° and 220° 
have been chosen that correspond to the "near" and "far" positions of the cross 
with respect to the camera. In principle, it is not necessary for the angle difference 
to be 180° and several rotation positions can be used instead of two. 
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Figure 13; The procedure for the source Z-position determination. 

For each of the projections, the position of the shadow crossing point along 
the camera axis has been calculated, as it is shown in Figurej^ The dependence 
of the shadow cross position on the phantom Z-position is presented in Figure p!4| 
The blue line corresponds to the measurements at 40° rotation angle and the red 
line corresponds to the measurements at 220° rotation angle. It is obvious that the 
crossing point of these two lines is the desired source position along the camera Zc 
axis in pixels. 
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Figure 14: Determination of the source position along the axis 


5.4 Determination of the source position along the Y^. axis 
and of the source to detector distance SDD 

To determine the source position along the Yc axis, the measurement scheme shown 


in Figure 15 can be used. One needs to perform several sets of measurements of 
the shadow coordinate obtained for different "cross" height positions by moving 
the camera towards and away from the source. Different "cross" height positions 
are achieved by simply rotating the phantom by various angles. 
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Camera position 



Figure 15: Determination of the source Yc position, and of the source to detector 
distance SDD 

In our case, the Yc coordinate of the shadow cross have been measured for 
three positions of the gantry rotation angle in five camera positions along the Xg 
axis connecting the camera to the source. In Figure the green line corresponds 
to the shadow position measurements at 30° rotation angle, blue line corresponds 
to 203° gantry rotation angle and the red line corresponds to 358° gantry rotation 
angle. All three lines have to cross at the same point, which is the position of the 
x-ray source along the camera Y^ axis. As it is seen from the graph, this method is 
suitable not only for determining of the source Y^. position, but also to determine 
the distance SDD from the source focal spot to the camera, which in this case 
turned out to be approximately two centimeters greater than the nominal. 
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Figure 16: Results of calculation of the source Yc position and correction of the 
SDD distance between the source and the camera 


6 The impact of the tomograph alignment to the 
image reconstruction quality 

Geometrical corrections obtained in the procedure of the tomograph alignment can 
be used as input parameters for the image reconstruction. The plastic phantom 
with two additional wires strung in its cross section has been used as a scanned 
sample. In the first case, the image reconstruction is performed using the nominal 
parameters of the geometry (Figure [T^. In the second case, corrections obtained 
in the procedure of the tomograph alignment have been applied (Figure [I^. 

From these figures it is evident that without geometrical corrections the 
reconstruction is of lower quality, with several artifacts. In particular, the wires 
on the reconstructed image are look duplicated and dashed. The implementation 
of the alignment corrections results in a significant improvement in the quality of 
the reconstructed image. 
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Figure 17: Phantom image reconstructed without corrections 



Figure 18; Phantom image reconstructed with the alignment corrections applied 


7 Spatial resolution 

The overall spatial resolution of the tomograph depends both on the quality of 
the original projections and also on the quality of the reconstruction. In the hrst 
case, it is influenced by the physical processes in the material of the detector, the 
detector pixel size, the size of the x-ray tube focal spot, and the zoom factor. Such 
factors as mechanical precision of the movement of the tomograph parts and the 
precise alignment of the system also affect the quality of the image reconstruc¬ 
tion. In this section several methods are described that were used to obtain a 


20 


quantitative estimation of the overall spatial resolution of the system. 

The hrst method is based on indirect estimation of the spatial resolution 
using a Gaussian approximation to the point spread function (PSF). Aluminum 
wire of the diameter about 2.7 mm has been used as the scanned object. To 
evaluate the spatial resolution of the tomograph, an image presenting the sum 
over 20 reconstructed slices of the wire is used (Figure [I^. A narrow area along 
the two axes presented by the red lines in the hgure, which can be considered as 
rectangular, is selected for htting. A function which is a convolution of a rectangle 
with a Gaussian is used here as an approximation function (Figure [20| . The 
rectangular function is a non-smeared initial image of the wire, while the Gaussian 
function is the system response to a point-like object. As a result of the ht, the 
value of (T = 34 pm for the Gaussian point spread function (PSF) is determined 
along both axes. (The reader may note that Figure 20 has some beam hardening 
visible as a cupping artefact, but this does not affect our results.) 
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Figure 19: Sum over 20 reconstructed slices of aluminum wire (the wire is slightly 
damaged from the one side). 
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Figure 20: Approximation of the wire cross-section using convolution of a Gaussian 
with a rectangular function along the Xg axis (top plot) and along the Yg axis 
(bottom plot). 

To get a direct estimate of the spatial resolution, three different samples of 
stainless mesh have been scanned (Figure [2l|. The hrst mesh has wires of 25 fxm 
diameter and 67 fim repeat distance, the second mesh has wires of 30 fim diameter 
and 80 /im repeat distance and the third mesh has wires of 50 /im diameter and 
265 fim repeat distance. Assuming that the response to a point object is Gaussian, 
the images of the two point-like objects become separated when the distance be¬ 
comes larger than 2a approximately. The smallest mesh is barely resolved, but for 
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the second sample (30 /im diameter and 80 /im repeat distance) the mesh struc¬ 
ture is fully visible. This is in a good agreement with the estimation above of the 
spatial resolution. 



Figure 21: Reconstructed images of the three samples of stainless mesh: 25 /rm 
diameter 67 /rm repeat distance, 30 /im diameter 80 fim repeat distance and 50 fim 
diameter 265 /rm repeat distance. 


A second independent direct check of a spatial resolution has been performed 
by scanning a standard calibration phantom, the QRM-Barpattern-Phantom (Fig¬ 


ure 


22). This calibration phantom consists of two perpendicular silicon 5x5 mm^ 


chips with pattern bars and holes of 5 to 150 /im size and of 80 to 120 /im depth 
(Figure [2^. Reconstructed images of the phantom are shown in Figure 23 The 
least size of visible bars on these pictures is 50 fim, this corresponds to 10 line-pairs 
per millimeter. 


Block 

linewidth (^Jm) 

linepairs per 
pattern 

points (pm) 

points per 
pattern 

A 

5, 10,25,50, 100, 150 

5 



B 

5, 10, 15, 20,25,30 

5 

5, 10, 15, 20, 25, 30 

18 

C 

5, 10, 15, 20,25,30 

5 

5, 10, 15, 20, 25, 30 

18 

D 



5, 10, 25, 50, 100, 150 

18 

E 

5,10,25,50, 100, 150 

5 




Figure 22: Description of the QRM-Barpattern-Phantom 
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Figure 23; Reconstructed images of QRM-Barpattern-Phantom. Left plot; recon¬ 
structed 3D image. Middle plot; phantom cross-section in the X-Y plane. Right 
plot; phantom cross-section in the X-Z plane. 


8 Conclusions 

In computer tomography, a sufficiently accurate knowledge of the tomograph ge¬ 
ometry is required to achieve good quality image reconstruction. In this paper, 
the precision of the movement of the MARS scanner parts has been studied and 
several geometrical parameters have been determined to improve the quality of 
the image reconstruction. Several procedures for measuring the MARS scanner 
alignment have been developed and applied. The resulting corrections were used 
as input parameters for the image reconstruction algorithm. The results of this 
work show that using these alignment corrections leads to a signihcant improve¬ 
ment in the quality of the image reconstruction. Using these corrections, several 
independent methods have been applied to make quantitative estimation of the 
hnal spatial resolution for the reconstructed image. The results of different meth¬ 
ods are consistent and give an estimation for the tomograph spatial resolution of 
slightly greater than 30 microns. 


References 

[1] MARS Bioimaing Ltd, 29a Clyde Road Christchurch, NZ 


25 





[2] L.Tlustosa, G.Shelkov, O. Tolbanov Characterisation of a GaAs(Cr) Medipix2 
hybrid pixel detector NIM Volume 633, Supplement 1, May 2011, Pages 
S103-S107 

[3] R. Ballabriga, M. Campbell, E. Heijne, X. Llopart, L. Tlustos, W. Wong 
Medipix3: : A 64 k pixel detector readout chip working in single photon 
counting mode with improved spectrometric performance Nucl. Instr. and 
Meth. A In Press, Corrected Proof, Available online 17 June 2010 

[4] RENISHAW RESR angle encoder system http;//www.renishaw.com 

[5] Langner 0.,Karolczak ,Rattmann and Kalender Bar and Point Test Patterns 
Generated by Dry-Etching for Measurement of High Spatial Resolution in 
Micro-CT; 2009 IFMBE Proceedings, World Congress on Medical Physics 
and Biomedical Engineering, September 7 - 12, 2009, Munich, Germany Vol. 
25/2, 428-431 Diagnostic Imaging 


26 



